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In this communication, wc consider the normal ordering of operators of the type 

a + = e+l 

where a (resp. a + )is a boson annihilation (resp. creation) operator; these satisfy 
[a, a+] = aa+ — a+a = 1, and for the purposes of this note may be thought of as 
a = d/dx and a~*~ = x. We discuss the integration of the onc-paramctcr groups 
e xn and their combinatorial by-products. In particular we show how these groups 
can be realized as groups of substitutions with prefunctions. 



1. Introduction 

This text is the continuation of a series of works on the combinatorial and 
analytic aspects of normal forms of boson strings^' 2 ' 3 ' 4,5 ' 11 ' 13 ' 14 ' 15 ' 18 ' 19 ' 20 ]. 
Let w G {a, a + }* be a word in the letters {a, a + } (i.e. a boson string), and 
define (as in Blasiak, Pcnson and Solomon 1 ' 2 ' 3 ' ) by r, s and e, respectively 
|to| Q + (the number of creation operators), |to| (the number of annihilation 
operators) and r — s (the excess) , then the normal form of w n is 



AT(w n ) = (a+r s ™( n > fc )(« + ) fc « fc (!) 

\fe=o / 

when e is positive (i.e. there are more creation than annihilation operators). 
In the opposite case (i.e. there are more annihilations than creations) the 
normal form of w n is 

JVV) = (f^S w (n 7 k)(a+) k a k J (a) n ^ (2) 
\fc=o / 

in each case, the coefficients S w arc defined by the corresponding equation 
(1 and 2). 

Now, for any boson string u one has 

Af(u) = (a+)l tt l £ .+ a |u| » + X " v - ( 3 ) 

\v\<\u\ 

It has been observed [ 12 ] that the numbers At, are rook numbers. 

Consider, as examples, the upper-left corner of the following (doubly 
infinite) matrices. 

For w — a + a, one gets the usual matrix of Stirling numbers of the 
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second kind. 



10 
10 



1 
1 
1 



00 

00 

00 

00 



1 



1 15 25 10 





1 



1 31 90 65 15 1 



For w = a + aa + , we have 



1 
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1 
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1 
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18 


9 


1 








24 


96 


72 


16 


1 





120 


600 


600 


200 


25 


1 



720 4320 5400 2400 450 36 1 



+ aaa + a + , 


one gets 












1 




















2 4 


1 

















12 60 


54 


14 


1 











144 1296 


2232 


1296 


306 


30 


1 






2880 40320 109440 105120 45000 9504 1016 52 1 



(4) 



(5) 



(6) 



Remark 1.1. In each case, the matrix S w is of staircase form and the 
"step" depends on the number of a's in the word w. More precisely, due to 
equation (3) one can prove that each row ends with a 'one' in the cell (n, nr), 
where r = \w\ a and we number the entries from (0, 0). Thus all the matrices 
are row- finite and unitriangular iff r = 1, which case will be of special inter- 
est in the following. Moreover, the first column is (1, 0, • • • , 0, • • • , 0, • • • ) 
iff w ends with an a (this means that Af(w n ) has no constant term for all 
n > 0). 



In this communication, we concentrate on boson strings and more gen- 
erally (homogeneous) boson operators involving only one "a" . We will see 
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that this case is closely related to one-parameter substitution groups and 
their conjugates. 

The structure of the paper is the following. 

In section 2 we define the framework for our transformation matrices 
(spaces, topology and decomposition), then we concentrate on the Riordan 
subgroup (i.e. transformations which are substitutions with prefunctions) 
and adapt the classical theory (Sheffer condition) to the present context. 
In section 3 we analyse unipotent transformations (Lie group structure and 
combinatorial examples). The divisibility property of the group of unipo- 
tent transformations tells us that every transformation is embedded in a 
one-paramater group. This will be analysed in section 4 from the formal 
and analytic points of view. Section 5 is devoted to some concluding re- 
marks and further interesting possibilities. 

2. The algebra £(C N ) of sequence transformations 

Let C N be the vector space of all complex sequences, endowed with the 
Frechet product topology. It is easy to check that the algebra £(C N ) of 
all continuous operators C N — ► C N is the space of row-finite matrices with 
complex coefficients. Such a matrix M is indexed by N x N and has the 
property that, for every fixed row index n, the sequence (M(n,k)) k>0 has 
finite support. For a sequence A — (a„)„> , the transformed sequence 
B = MA is given by B = (b n ) n >o with 



Remark that the combinatorial coefficients S w defined above are indeed 
row-finite matrices. 

We may associate a univariate series with a given sequence (o n )„gNi us- 
ing a sequence of prescribed (non-zero) denominators (cZ„) rae N ; as follows: 



For example, with d n = 1, we get the ordinary generating functions (OGF), 
with d n = n\, we get the exponential generating functions (EGF) and with 
d n = (nl) 2 , the doubly exponential generating functions (DEGF) and so on. 
Thus, once the denominators have been chosen, to every (linear continuous) 
transformation of generating functions, one can associate a corresponding 
matrix. 




(7) 



fe>0 




(8) 
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The algebra £(C N ) possesses many interesting subalgebras and groups, 
such as the algebra of lower triangular transformations T(N, C), the group 
T inv (N, C) of invertible elements of the latter (which is the set of infinite 
lower triangular matrices with non-zero elements on the diagonal), the sub- 
group of unipotent transformations UT(N, C) (i.e. the set of infinite lower 
triangular matrices with elements on the diagonal all equal to 1) and its Lie 
algebra WT(N, C), the algebra of locally nilpotcnt transformations (with 
zeroes on the diagonal). One has the inclusions (with 2?, n „(N, C), the set 
of invertible diagonal matrices) 



We remark that T mv (N, C) = X> in „(N, C) cx UT(N, C) because UT is 
normalized by V inv and % nv — V inv .UT (every invertible transformation 
is the product of its diagonal by a unipotent transformation). 

We now examine an important class of transformations of T as well as 
their diagonals: the substitutions with prefunctions. 

2.1. Substitutions with prefunctions 

Let (d n ) n >o bet a fixed set of denominators. We consider, for a generating 
function /, the transformation 



The matrix of this transformation M g ^ is given by the transforms of the 
monomials §- hence 



UT(N, C) c T inv (N, C) c T(N, C) c £(C N ) 
V inv (N,C) c 7^(N,C) and WT(N, C) c £(C N ). 



(9) 



QgAfK*) = atom*))- 



(10) 




(11) 



If 



g, 4> ^ (otherwise the transformation is trivial), we can write 




(12) 



with ai 7 a m ^ and then, by (11,12) 




(13) 



One then has 
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M g ^ is row — finite <==^> <p has no constant term (14) 
and in this case it is always lower triangular. 

The converse is true in the following sense. Let T e £(C N ) be a matrix 
with non-zero two first columns and suppose that the first index n such 
that T(n, k) ^ is less for k = than k = 1 (which is, from (11) the case 
when T = M g ^). Set 

3 (x):=d £T(n,0)^; <j>{x) := £ T(n, 1)- (15) 

n>0 ™ 5 ^ ' n>0 ™ 

then T = Mgj iff, for all k, 

n>0 

Remark 2.1. Eq. (11) is called the Sheffer condition (see [i6,20,2i,23j^ anc ^ 
for EGF (d n = n\) it amounts to stating that 



V T(n,k)^-y k =g(x)e^ x \ (17) 
' n! 



Tl.fcX) 



From now on, we will suppose that <ft has no constant term (ao = 0). 
Moreover M g ^ <G 7^ if and only if a , a.\ ^ and then the diagonal term 

with address (n,n) is |£ . We get 

M,^WT^^ = ^ = 1. (18) 

In particular for the EGF and the OGF, we have the condition that 

g(x) = 1 + higher terms and <j>(x) = x + higher terms. (19) 

Note 2.1. In classical combinatorics (for OGF and EGF), the matrices 
M g rf{n,k) are known as Riordan matrices (see [ 16 » 17 ] for example). 

3. Unipotent transformations 
3.1. Lie group structure 

We first remark that n x n truncations (i.e. taking the [0..n] x [0..n] sub- 
matrix of a matrix) are algebra morphisms 

r„ : T(N, C) -» M([0..n] x [0..n], C). (20) 
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We can endow T(N, C) with the Frechet topology given by these mor- 
phisms. We will not develop this point in detail here, but this topology is 
metrisable and given by the following convergence criterion: 

a sequence (M&) of matrices in T(N, C) converges iff 

for all fixed n G N 
the sequence of truncated matrices (r„(Mfe)) converges. 

This topology is compatible with the C-algebra structure of T(N, C). 

The two maps exp : WT(N, C) -» «T(N, C) and log : WT(N, C) -» 
A/ r T(N, C) are continous and mutually inverse. 

3.2. Examples 

3.2.1. Provided by the exponential formula 

The "classical exponential formula" [ 7 > 9 < 21 ] tells us the following: Consider 
a class a of finite labelled graphs C. Denote by C c the subclass of connected 
graphs in C. Then the exponential generating series of C and C c are related 
as follows: 



The following examples give us some insight into why combinatorial matri- 
ces of the type: 

T(n, k) = Number of graphs of C on n vertices having k connected components 
give rise to substitution transformations. 
Example 3.1. Stirling numbers. 

Consider the class of graphs of equivalence relations. Then using the statis- 
tics ^( num ^ cr °* points) ^.(number of connected components) 



EGF(C) = e EGF(cc K 



(21) 



E S{n,k)-y 



n,k>0 




X 



(number of points of T) 



(number of connected components of T) _ 



(number of points of T)\ 



y 



all equivalence graphs T 



a Closed under relabelling (of the vertices), disjoint union, and taking connected compo- 
nents. 
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^(number of vertices of T) 



exp 



^- — 4 (number of points of T)! 

*X connected 



-y 



(number of connected components of F) 



exp 



E^ 

,n>l 



= e y(e x -i) 



(22) 



we will see that the transformation associated with the matrix S(n, k) 
is / - f(e* - 1). 

Example 3.2. Idempotent numbers. 

We consider the graphs of endofunctions b . Then, using the statis- 
tics ^( num ^ cr °f points of the set)^(numbcr of connected components of the graph) 

denoting by I(n, k) the number of endofunctions of a given set with n 
elements having k connected components, we get 



E I{n,kf-y* = 



n,fc>0 



E 



all graphs of endofunctions T 

exp 



^.(number of vertices of T) 
(number of vertices of T)l 



(number of connected components of T) 



r (number of vertices of T) 

(number of connected components of T) 

, y 



\T connected 



(number of vertices of Y)\ 



exp 



E^ 



, n>l 



(23) 



Corresponding to these numbers we get the (doubly) infinite matrix 



10 

1 

02 1 00 

03 6 1 00 
4 24 12 1 
5 80 90 20 1 
6 240 540 240 30 1 



(24) 



b Functions from a finite set into itself. 
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and wc will see that the transformation associated with this matrix is 



3.2.2. Normal ordering powers of boson strings 

To get unipotent matrices, one has to consider boson strings with only one 
annihilation operator. In the introduction, we have given examples with 
a + a,a + aa + (the matrix of the third string, a + aaa + a + , with two annihila- 
tors, is not unipotent). Such a string is then w = (a + ) r ~ p a(a + ) p and we 
will see shortly that 

• if p = 0, S w (n, k) is the matrix of a unipotent substitution 

• if p > 0, S w (n,k) is the matrix of a unipotent substitution with 
prefunction 

To cope with the matrices coming from the normal ordering of powers 
of boson strings we have to make a small detour to analysis and formal 
groups. 

4. One-parameter subgroups of UT(N, C) 
4.1. Exponential of elements of NT(N, C) 

Let M = I + N e UT(N, C) (7 = J N is the indentity matrix). One has 



One can see that, for k < n, due to the local nilpotency of N, the ma- 
trix coefficient M'(n, k) is well defined and, in fact, a polynomial of degree 
n — k in t (for k > n, this coefficient is 0). We have the additive property 
]\/fti+t 2 _ M tr M t2 and the correspondence t — » M* is continuous. Con- 
versely, let t — > M t be a continous local one-parameter group in UT(N, C); 
that is, for some real e > 



/- fixe") 




(25) 



where 





(26) 



|*i | and \t 2 \ < e => M tl M ±2 = M tl +t 2 



(27) 
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then there exists a unique matrix H e NT(N, C) such that M t — 
exp(tH). (This may be proved using the projections t„ and the classical 
theorem about continuous one-parameter subgroups of Lie groups, see [ 10 ], 
for example). 

When M t — M* is defined by formula (25) we have 
H = log(I + N) = j: k > 1 ^f 1 N k . 

The mapping t — > M* will be called a one parameter group of UT(N, C). 



Proposition 4.1. Lei M 6e i/ie matrix of a substitution with prefunction; 
then so is M* for all t e C. 

The proof will be given in a forthcoming paper and uses the fact that 
"to be the matrix of a substitution with prefunction" is a property of poly- 
nomial type. But, using composition, it is straightforward to see that M* 
is the matrix of a substitution with prefunction for all teN. Thus, using 
a "Zariski-type" argument, we get the result that the property is true for 
all t e C. 



4.2. Link with local Lie groups : Straightening vector fields 
on the line 

We first treat the case p = of subsection (3.2.2). The string (a + ) r a cor- 
responds, in the Bargmann-Fock representation, to the vector field x r 
defined on the whole line. 

Now, we can try (at least locally) to straighten this vector field by a dif- 
feomorphism u to get the constant vector field (this procedure has been 
introduced by G. Goldin in the context of current algebras[ 8 ]). As the one- 
parameter group generated by a constant field is a shift, the one-parameter 
(local) group of transformations will be, on a suitable domain 

U x [f](x) = f(u- 1 (u(x) + X)). (28) 

Now, we know from section (4.1) that, if two one-parameter groups have the 
same tangent vector at the origin, then they coincide (tangent paradigm). 
Direct computation gives this tangent vector : 



d 

dX 



f(u-'(u(x) + X)) = ^ r -f(x) (29) 

A=0 U \ X ) 



and so the local one-parameter group U\ has -^rj^-£: as tangent vector 
field. 
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Here, we have to solve u , 1 ^ = x r in order to get the diffeomorphism u. 
In the case r ^ 1, we have (with V =]0, +oo[ as domain) 

1— r 

u(x) = ^—=y:u- 1 (y) = ((l-r)y)^ (30) 

and 

e Xxr Hf]{x) = f ( — —A (31) 

A(r - l)^- 1 )— / 

The substitution factor s\(x) = - r— has been already 

obtained by other means in The computation is similar for the case 
when r = 1 and, for this case, we get 

e^%](z) = /(eV) (32) 

with s\(x) — e x x as substitution factor. 

These first examples are summarized in the following table 



r 


s\(x) 


Name 





x + A 


Shift 


1 


e x x 


Dilation 


2 


X 

l-\x 


Homography 


3 


X 

Vl-2\x 2 





Comment 4.1. If one uses classical analysis (i.e. convergent Taylor series), 

one must be careful about the domain where the substitutions arc defined 

and the one-parameter groups are defined only locally. 

For each of these examples, one can check by hand that for suitable (and 

small) values of A, one has s\{s^(x)) = s\+^{x) (one-parameter group 

property). 

It is possible to avoid discussion of the domains by considering A, /i as new 
variables and applying the "substitution principle"; namely by claiming 
that it is possible to substitute a series without constant term in a series 
(in the algebra C[[x, A, fj]]). 

Using the same method, one can start with more complicated operators. 
Examples and substitution factors are given below 
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Operator 


Substitution Factor 


Description 


(l + (a+) 2 )a 


( ^ xcos(X) + sin(X) 
SX X cos(X) — xsin(X) 


One-parameter group 
of homographies 


v /l + (a+) 2 
a+ 






s x (x) = \/x 2 + 2AV1 + x 2 + X 2 


Composition of quadratic 




direct and inverse functions 



4.3. Case p > 0; another conjugacy trick and a surprising 
formula 

Regarding vector fields as infinitesimal generators of one-parameter groups 
leads to conjugacy since, if U\ is a one-parameter group of transformation, 
so too is VUxV^ 1 (V being a continuous invertible operator). We could 
formally consider (a + ) r ~ p a(a + ) p ; p > as conjugate to ((a + ) r a) in this 
context. More generally supposing all the terms well-defined, if 

fi = ui(x)—u 2 (x) = — — ( ui(x)u 2 (x)—) u 2 (x) 
ax u 2 (X) \ ax J 



then 



e \n = 1 f\ Ul {x)u 2 (x)± 
u 2 {x) V 



) u 2 (x) (33) 



This rather surprising formula (33) may be understood as an operator 
equality. 

Now, the tangent paradigm (see section 4.2) tells us that, if we adjust this 
tangent vector to coincide with x r ~ p -^x p (recall that the original problem 
was the integration of the operator Q = (a + ) r ~ p a(a + ) p ; p > 0), then we 
get the right one-parameter group. Using this "conjugacy trick" we get 

e xn [f](x) = ( S -^-) f(s x (x)) with s x (x) = -3- (34) 



(1 - X(r - l)x r 

Remark 4.1. (i) ft can be checked that, if s\{x) is a substitution factor 
(i.e. at least locally s\(s l _ L (x)) = s\ + ^(x)) such that s\(0) = for every 
A (which is the case in most of our examples) then the transformations 
defined by U\[f](x) = ^ Sx ^ f(s\(x)) form a one-parameter (possibly 
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local) group. 

(ii) It is also possible to use the "ad" operator (Lie adjoint) instead of "Ad" 
(conjugacy) to obtain integration formulas (see Dattoli 6 ). 

4.4. Characteristic series «-> one parameter group 
correspondence 

The preceding allows us to extend integration processes to linear combina- 
tions of boson strings in the following sense. The algebra W\ j(yo generated 
by a + , (a + ) _1 , a is graded by 

weight(a+) = 1, weight ((a + ) _1 ) = weight(a) = -1 (35) 

and every homogeneous operator of this algebra which is of the form 

O = a w w (36) 

|w|a = 1 5 wcight(^)— e 

(there is only one annihilation operator in each monomial) can be integrated 
as above. So one would like to reconstruct the characteristic series 

J2Mn,k)—y k (37) 

n,k 

from knowledge of the one-parameter subgroup e xn . 
This is the aim of the following paragraph. 

For every homogeneous operator as above with e > 0, one defines the 
coefficients Sn(n, k) as in the Introduction by 

oo 

W(Q") = (a+) ne Mn, k){a+) k a k (38) 

fe=0 

One has the following proposition 

Proposition 4.2. With the definitions introduced, the following conditions 
are equivalent: 

V Sn(n,k)^-y k =g(x)eyt>W (39) 
z — ' n! 

n,k>0 

Ux[f](x) = g(Xx e )f (x (1 + <HAx e ))) (40) 
Which solves the problem. 
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5. Conclusion and remaining problems 

We have considered a class of elements of Wi )00 (see section 4.4 for a def- 
inition) which describe some rational vector fields on the line. For these 
operators, we have established a correspondence 

One-parameter group (integration of the field) <-> 
Characteristic series (^coefficients of the normal ordering) 

We have then seen that families of combinatorial matrices give rise through 
the exponential formula to substitutions. 

Further work which remains is to study the vector fields associated with 
these combinatorial matrices. Also it would be desirable to adapt this 
machinery to other algebras (quons, several boson modes). 
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